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ABSTRACT 


Autonomous  vehicles  (AVs)  are  commonly  used  in  oceanic  and  more  recently  estuarine 
and  riverine  environments  because  they  are  small,  versatile,  efficient,  moving  platforms 
equipped  with  a  suite  of  instruments  for  measuring  environmental  conditions.  However, 
moving  vessel  observations,  particularly  those  associated  with  Acoustic  Doppler  Current 
Profiler  (ADCP)  measurements,  can  be  problematic  owing  to  instrument  noise,  flow 
fluctuations,  and  spatial  variability.  A  range  of  ADCPs  manufactured  by  different 
companies  were  integrated  onto  an  Unmanned  Surface  Vehicle  (USV),  an  Unmanned 
Underwater  Vehicle  (UUV),  and  some  additional  stationary  platforms,  and  were 
deployed  in  a  number  of  natural  riverine  and  estuarine  environments  to  evaluate  the 
quality  of  the  velocity  profile  over  the  depth,  minimum  averaging  time  interval 
requirements  and  AV  mission  planning  considerations.  An  appropriate  averaging 
window,  T*,  was  detennined  using  the  Kalman  Algorithm  with  a  Kalman  gain  equal  to 
1%.  T*  was  found  to  be  independent  of  depth,  flow  velocity,  and  environment.  There  was 
no  correlation  (R  =0.18)  for  T*  between  flow  magnitude  and  direction.  Results  from  all 
measurements  had  a  similar  T*  of  approximately  3  minutes.  Based  on  this,  an  averaging 
window  of  four  minutes  is  conservatively  suggested  to  obtain  a  statistically  confident 
measure  of  the  mean  velocity  profile. 
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I.  INTRODUCTION 


The  versatility  and  flexibility  of  autonomous  vehicles  (A Vs)  provide  a  unique 
environmental  survey  platform.  AVs  are  small  in  size,  typically  capable  of  being 
deployed  and  operated  by  one  person,  and  are  equipped  with  a  sensor  suite  comparable  to 
those  generally  mounted  on  larger  sized  vessels.  The  size  and  weight  of  sensors  continue 
to  decrease  while  vehicle  functionality  and  capability  continue  to  increase,  providing 
scientists  with  a  new  set  of  tools  to  improve  our  understanding  of  various  aquatic 
environments.  In  addition  to  shrinking  the  size  of  vehicle  platforms  and  sensors,  the 
associated  cost  of  AVs  are  becoming  more  affordable,  especially  when  the  operating 
costs,  such  as  vessel  maintenance,  personnel,  and  fuel,  are  included.  This  affords  the 
scientist  with  the  ability  to  procure  and  deploy  AVs  in  natural  environments,  allowing  for 
greater  spatial  coverage,  reducing  the  time  necessary  to  complete  data  collection,  and 
reducing  logistical  costs.  AVs  come  in  a  variety  of  shapes  and  sizes  resulting  in  different 
capabilities  and  limitations. 

Standard  aquatic  AVs  are  equipped  with  a  combination  of  positioning,  depth, 
velocity,  and  water  quality  sensors.  Since  most  AVs  are  continuously  moving,  issues  of 
instrument  noise,  environmental  noise  (scales),  and  stationarity  (if  averaging  is 
perfonned)  must  be  considered.  Instrument  noise  is  assumed  to  be  a  random  process. 
Therefore,  statistical  confidence  of  the  mean  is  gained  by  obtaining  a  number  of 
observations,  which  are  then  averaged.  The  greater  the  number  of  independent 
observations  ( N)  collected,  the  higher  the  confidence  in  the  mean  (Bendat  and  Piersol 
2000).  However,  there  is  a  point  of  diminished  returns,  and  detennining  this  point  (in 
time)  is  critical  for  designing  an  efficient  sampling  scheme.  The  number  of  observational 
points  considered  to  be  independent  is  given  by 


N 


T 


t 


decorr 


n-  dt 


v  decorr 


(1) 


where  tdecon-  is  the  environmental  decorrelation  time,  T  is  the  total  sampling  duration,  and 
T  =  n-  dt,  where  n  is  the  total  number  of  observations  and  dt  is  the  instrument  sampling 
rate.  For  random  data  (e.g.  instrument  noise)  tdecon ■  =  dt,  and  each  observational  data  point 
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is  considered  independent  regardless  of  sampling  frequency  ( N=n ).  However,  in  nature 
there  can  be  random  processes  that  have  statistically  ensemble-defined  temporal  and 
spatial  scales  leading  to  longer  decorrelation  times.  For  example,  if  tdecorr=  10  s,  and  an 
instrument  samples  for  7=  1 00  s,  the  record  will  consist  of  N=  1 0  independent 
observations,  or  degrees  of  freedom,  regardless  of  the  instrument  sample  rate.  If  the 
observations  are  collected  from  a  moving  platform,  spatial  limitations  are  also  required 
such  that  the  mean  and  environmental  scales  (noise)  are  not  evolving  (stationarity)  within 
the  sampling  duration  required  to  obtain  a  statistically  confident  estimate  of  the  mean. 
Therefore,  the  appropriate  AV  sampling  scheme  is  dependent  upon  instrument  noise,  the 
environment,  the  vessel  speed,  and  the  observations  of  interest  including  the  required 
statistical  confidence.  In  order  to  ensure  a  true  estimate  of  the  mean  velocity  is  obtained 
when  evaluating  an  appropriate  averaging  window,  and  to  eliminate  errors  due  to  spatial 
variability  from  moving-vessel  measurements,  only  measurements  obtained  by  the 
vessels  at  fixed  locations  are  considered  in  this  work. 

Acoustic  Doppler  Current  Profilers  (ADCPs)  were  originally  implemented  by 
scientists  on  aquatic  AVs  to  provide  dead-reckoning  navigation  based  on  a  Doppler 
Velocity  Log  (DVL)  (Fong  and  Jones  2006).  However,  ADCPs  have  recently  been  used 
on  AVs  to  measure  features  including  coastal  ocean  flows  (An  et  al.  2001;  Dhanak  et  al. 
2001;  Fong  and  Jones  2006;  Hibler  et  al.  2008;  Shay  and  Cook  2003),  directional  surface 
waves  (Hayes  et  al.  2007)  and  depth-averaged  currents  (Eriksen  et  al.  2003).  There  is  an 
increasing  need  for  collecting  environmental  data  with  AVs  in  faster  and  more  dynamic 
flows  found  in  riverine  and  estuarine  environments.  Historically,  in  rivers,  ADCPs  have 
been  mounted  on  moving  vessels  to  measure  discharges  (Gordon  1989;  Morlock  1996; 
Muste  et  al.  2004a;  Yorke  and  Oberg  2002),  but  there  is  a  growing  interest  in  collecting 
ADCP  measurements  for  describing  mean  velocity  profiles  over  the  depth.  Sophisticated 
data  analysis  techniques  have  been  developed  that  are  able  to  resolve  the  depth-averaged 
tidal  and  non-tidal  flow  behavior  from  moving-vessel  (noisy)  ADCP  measurements 
(Vennell  and  Beaston  2006;  Vennell  and  Beaston  2009),  however  the  focus  here  is  to 
resolve  depth-varying  currents.  Current  profiles  have  been  successfully  measured  using 
fixed  ADCPs  in  oceanic  and  estuarine  environments,  where  ADCPs  were  either 
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downward  looking  from  buoy  moorings  or  from  the  hulls  of  moored  ships,  or  upward 
looking  from  bottom  mounts  (Gordon  et  al.  1990). 

ADCP  measurements  are  inherently  noisy.  Muste  et  al.  2004b  recommend  that 
fixed  vessels  collect  7  to  1 1  minutes  of  stationary  data  in  order  to  resolve  the  instrument 
and  environmental  noise  associated  with  ADCP  measurements  in  riverine  environments. 
However,  this  recommended  averaging  time  is  based  on  visual  estimates  of  when  the 
mean  flow  velocity  becomes  stable  using  data  collected  in  only  one  environment,  and 
therefore  requires  further  examination.  In  order  to  accurately  capture  mean  flow 
characteristics,  the  relationship  between  ADCP  flow  measurements  obtained  with  A  Vs  in 
differing  non-wavy  environments  and  the  ADCP  operational  parameters  also  need 
additional  evaluation.  This  paper  focuses  on  two  different  AVs  (SeaRobotics  USV-2600 
Unmanned  Surface  Vehicle,  and  YSI/Oceanserver  Iver2  Unmanned  Underwater  Vehicle) 
equipped  with  ADCPs  that  were  deployed  in  a  river  and  tidal  inlet  and  collected 
measurements  at  fixed  locations,  with  the  primary  objective  of  obtaining  statistically 
confident  depth-varying  velocity  profiles  and  detennining  the  associated  averaging 
window  necessary  to  remove  instrument  and  environmental  noise,  as  discussed  by  Muste 
et  al.  (2004b).  Knowledge  of  the  appropriate  averaging  window  required  to  obtain  a 
statistically  confident  measure  of  the  mean  velocity,  determined  by  robust  statistical 
techniques,  allows  users  to  maximize  their  AV  mission  planning  to  optimize  time  and 
spatial  resolution.  Proper  mission  planning  is  an  essential  part  of  protecting  U.S.  Navy 
personnel  and  equipment  from  unnecessary  risks.  This  type  of  data  collection  requires  the 
users  and  equipment  to  remain  stationary  for  relatively  long  periods  of  time.  Therefore, 
Navy  personnel  and  equipment  are  subject  to  periods  of  exposure  to  any  threats  that 
particular  environment  may  possess.  Additionally,  proper  mission  planning  will  provide 
the  ability  for  data  collection  over  larger  areas,  which  provides  the  Navy  with  a  greater 
understanding  of  the  environment  and  increased  confidence  for  environmental 
characterization.  Thus,  determining  the  appropriate  mission  parameters  is  considered  a 
critical  need  for  the  U.S.  Navy. 
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II.  AUTONOMOUS  VEHICLES 


A.  UNMANNED  SURFACE  VEHICLE  (USV) 

The  Sea  Robotics  USV-2600  (Figure  1),  discussed  herein,  is  a  1.9  x  1.25  m 
fiberglass  twin-hull  catamaran,  which  weighs  about  360  lbs  with  all  instrumentation 
included.  It  is  powered  by  dual  weedless  propellers,  located  about  0.5  m  below  the  water 
surface.  The  USV  is  designed  to  support  multiple  in-situ  sensors  that  can  be  configured  to 
meet  the  user’s  requirements  for  taking  measurements  in  any  type  of  environment.  The 
USV  supports  a  1.2  MHz,  4-beam  bottom-tracking  RD Instruments  (RDI)  ADCP,  a 
Differential  Global  Position  System  (DGPS)  antenna,  a  single-beam  echosounder, 
heading,  altitude,  and  several  water  quality  sensors.  It  has  an  onboard  data  storage 
capacity  of  3  GB  and  provides  the  user  plenty  of  memory  to  complete  a  full  day  of  data 
collection.  Data  can  also  be  streamed  real-time  via  radio  antenna  over  the  2.4  GHz 
frequency  band  as  long  as  line-of-site  with  the  USV  is  achieved.  The  USV  is  powered  by 
two  Lithium  Polymer  battery  packs  capable  of  producing  a  top  speed  of  4.7  m/s.  The 
USV  can  survey  a  given  environment  for  up  to  10  hours  while  operating  at  1.5  m/s  and 
can  last  up  to  1.5  hours  at  3  m/s. 

The  SeaRobotics  USV-2600  has  the  unique  capability  to  station-keep  owing  to  its 
dual-propeller  design.  The  station-keeping  feature  is  the  best  method  for  collecting 
statistically  reliable  data  as  it  removes  the  spatial  sampling  constraints,  e.g.  stationarity. 
The  USV  does  not  have  an  Inertial  Measurement  Unit  (IMU),  however  since  USVs 
operate  on  the  water  surface,  continuous  DGPS  allow  for  precise  navigation  and  position 
recording,  and  vehicle  motion  is  accounted  for  in  the  ADCP  data  processing  associated 
with  bottom-tracking.  USVs  tend  to  be  larger  than  their  Unmanned  Underwater  Vehicle 
counterparts,  requiring  a  wheeled-trailer  for  transport  and  deployment.  The  USVs 
increased  size  tends  to  provide  an  increase  in  payload,  allowing  for  more  battery  capacity 
and  increasing  operational  speed  and  duration.  Making  additions  and  modifications  to 
USVs  also  tend  to  be  easier  for  the  end-user  owing  to  the  large  deck  space,  which 
provides  flexibility  to  adapt  the  USV  for  various  environmental  measurement  needs. 
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B.  UNMANNED  UNDERWATER  VEHICLE  (UUV) 

The  YSI/Oceanserver  EcoMapper  Iver2  UUV  (Figure  2),  discussed  herein, 
measures  1.6  m  long  with  a  diameter  of  0.15  m  and  weighs  45  lbs  in  air.  The  UUV  can 
operate  at  depths  down  to  60  m  using  four  independent  control  planes,  and  can  travel  at  a 
speed  of  0.5  to  2  m/s,  with  a  maximum  speed  of  1  m/s  on  the  surface.  For  navigation  at 
the  surface  and  to  verify  its  position,  the  UUV  uses  GPS  with  Wide  Area  Augmentation 
System  (WAAS)  corrections  to  provide  positional  accuracy  better  than  3  meters.  The 
UUV  configuration  includes  a  Sontek  10-Beam,  bottom-tracking  DVL  (up-  and  down¬ 
looking  configurations  consisting  of  four  velocity  beams  operating  at  1.0  MHz  and  a 
vertical  center  beam  operating  at  0.5  MHz  each),  allowing  for  improved  underwater 
navigation  and  water  current  profiling,  a  dual-frequency  side-scan  sonar  for  bottom 
imaging,  and  a  full  suite  of  water  monitoring  sensors  that  measure  conductivity, 
temperature,  depth,  pH,  dissolved  oxygen,  chlorophyll,  blue-green  algae,  turbidity,  and 
rhodamine.  The  UUV  does  not  have  an  IMU,  however,  the  vehicle  motion  is  accounted 
for  in  the  ADCP  data  processing  associated  with  bottom-tracking.  The  UUV  has  an 
onboard  data  storage  capacity  of  10  GB  and  runs  on  rechargeable  lithium-ion  batteries, 
making  it  capable  of  up  to  8  hours  of  data  collection  at  a  speed  of  1.3  m/s  in  a  no  flow 
environment. 

UUVs  are  now  stable  and  reliable  platforms  for  conducting  continuous  surveys  of 
the  water  column  in  environments  with  flows  less  than  a  few  knots.  Unlike  US  Vs,  UUVs 
are  normally  equipped  with  a  single  propeller,  which  limits  the  UUV’s  turning  radius  and 
its  ability  to  station-keep.  The  YSFOceanserver  EcoMapper  Iver2  UUV  does  have  the 
ability  to  remain  quasi-stationary  by  using  a  “park”  mode  for  a  user-defined  duration.  The 
“park”  mode  is  only  functional  at  the  surface  and  allows  the  UUV  to  station-keep  at  a 
specific  point  within  a  user-defined  radius  utilizing  its  GPS  antenna  with  WAAS 
corrections.  The  UUV  “park”  mode  is  defined  by  an  inner  and  outer  park  radius,  with  a 
minimum  outer  radius  of  3  m.  In  a  unidirectional  current,  the  vehicle  drives  to  the  center 
of  the  inner  park  radius  then  turns  off  the  propeller  and  floats.  When  the  vehicle  floats 
downstream  beyond  the  outer  park  radius,  it  drives  back  to  the  inner  park  radius, 
repeating  this  for  the  defined  park  duration.  Field  experience  has  shown  that  the  “park” 
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mode  works  best  when  the  vehicle  is  traveling  upstream  and  in  moderate  to  fast  flows 
and  can  hold  position  within  a  10  m  linear  excursion.  Although  outside  of  the  original 
intent  of  the  park  mode,  we  utilized  this  mode  for  obtaining  a  statistically  confident 
velocity  profile  at  a  relatively  stationary  location. 

There  are  several  logistical  problems  with  using  the  UUV  in  park  mode  because  it 
must  operate  at  the  surface.  Since  the  UUV  is  small  and  a  majority  of  the  body  floats 
below  the  surface,  it  can  be  difficult  to  see  by  other  boat  operators.  For  this  reason,  it  is 
recommended  that  a  support  vessel  be  used  to  warn  boat  traffic  about  the  UUV’s 
location.  A  tethered  surface  float  can  be  attached  to  the  UUV  to  increase  visibility,  but 
this  increases  the  drag  on  the  UUV  and  decreases  operational  time.  Also,  the  effect  of 
short  period  wind  waves  may  cause  the  UUV  to  pitch  and  roll  significantly  more  while  at 
the  surface.  Lastly,  but  most  importantly,  the  probability  of  floating  surface  debris,  such 
as  grass  and  seaweed,  getting  wrapped  around  the  propeller  increases  when  the  UUV 
operates  at  the  surface.  It  is  recommended  that  the  UUV  operate  below  the  surface  at  >1 
m/s  where  the  vehicle  is  stable  and  less  prone  to  propeller  fouling  and  being  hit  by  other 
boaters.  However,  caution  is  required  when  analyzing  the  fast-moving  vessel  ADCP  data 
to  ensure  environmental  flow  characteristics  are  not  changing  over  the  spatial  range 
covered  during  the  averaging  duration. 
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Figure  1.  SeaRobotics  USV-2600. 


Figure  2.  YSFOceanserver  EcoMapper  Iver2  UUV. 
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III.  EXPERIMENTS 


The  A  Vs  equipped  with  ADCPs,  as  well  as  several  other  platforms  with  velocity 
profiling  instruments,  were  deployed  in  various  environments  under  differing  flow 
conditions  (Table  1)  to  evaluate  the  performance  of  the  AVs  and  to  detennine  the 
appropriate  sampling  techniques. 

A.  KOOTENAI  RIVER,  ID,  AUGUST  2010 

ADCP  velocity  profiles  were  collected  by  the  USV  as  part  of  a  riverine  field 
experiment  conducted  in  August  2010,  on  the  Kootenai  River,  ID,  referred  to  as  KR 
(Figure  3A-C).  The  primary  goal  of  the  experiment  was  to  accurately  measure  the  3D 
flow  field  in  a  natural  river  composed  of  varying  depths  and  channel  meanders.  The 
backwater  meandering  reach  of  the  KR  (Figure  3B),  which  was  approximately  8  m  deep, 
200  m  wide,  and  had  flows  of  0.4  m/s,  was  measured  by  the  USV.  The  meandering  reach 
was  divided  into  14  transects  oriented  nonnal  to  the  river  bank  with  variable  streamwise 
spacing  for  enhanced  resolution  of  the  flow  dynamics  around  the  river  bends.  Each 
transect  consisted  of  five  locations  spaced  equally  across  the  river,  at  which  the  USV 
would  station-keep  to  within  a  few  meters  for  approximately  10  minutes.  The  ADCP 
mounted  on  the  USV  sampled  at  0.5  Hz,  and  acquired  velocities  throughout  the  water 
column  with  a  surface  blanking  distance  of  0.35  m,  0.25  m  bin  size,  and  28  depth  bins. 
Precise  positioning  and  navigation  to  each  transect  and  subsequent  stationary  profile 
location  were  achieved  by  the  onboard  DGPS.  An  additional  onboard  survey-grade 
DGPS  was  post-processed  after  the  mission  to  improve  the  USV  positions.  It  was  found 
that  the  USV  could  maintain  positioning  to  within  1  m.  This  high  degree  of  positional 
accuracy  is  necessary  for  describing  the  complex  flow  structures  that  can  occur  across  the 
width  of  a  river.  The  efficiency  of  the  USV  to  go  to  a  specified  location,  station-keep,  and 
then  continue  cannot  be  matched  by  single-propeller,  human-controlled  vessels,  which 
generally  have  the  additional  requirement  of  anchoring. 

In  addition,  self-contained,  downward-looking  2  MHz  Nortek  Aquadopp  ADCPs 
mounted  on  surface,  non-motorized,  mini-catamarans  (Figure  4)  were  deployed  in  the 
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braided  reach  of  the  KR  (Figure  3C),  which  was  3  m  deep,  100  m  wide,  and  had  flows  of 
1.5  m/s.  Three  mini-catamaran- ADCP  systems,  equally  spaced  in  the  across-stream 
direction,  were  hand  deployed  from  a  boat  in  transects  along  the  braided  reach.  At  each 
transect  the  mini-catamarans  were  anchored  to  the  bottom  for  10  minute  durations  to 
collect  stationary  ADCP  measurements.  The  ADCPs  on  the  mini-catamarans  sampled  at 
1  Hz,  and  had  a  surface  blanking  distance  of  0.05  m,  0.2  m  bin  size,  with  35  depth  bins. 

B.  ELKHORN  SLOUGH,  CA,  NOVEMBER  2010 

A  subsequent  6-hr  USV  mission  was  performed  in  November  2010,  in  Elkhorn 
Slough,  Monterey  Bay,  CA,  referred  to  as  ES,  which  is  a  shallow,  tidally-driven  slough 
(Figures  3D  and  3E).  ES  is  approximately  10  km  long  and  consists  of  a  main  channel 
with  a  complex  curving  structure,  mud  flats,  a  salt  marsh  and  numerous  small  tidal 
channels  (Breaker  et  al.  2003).  The  ADCP  on  the  USV  sampled  at  0.5  Hz,  and  measured 
velocities  over  the  water  column  with  a  surface  blanking  distance  of  0.35m,  a  0.25  m  bin 
size,  and  40  depth  bins.  The  USV  collected  data  over  ebb  and  flood  tidal  conditions, 
station-keeping  six  times  at  the  same  point  for  30  minutes  at  a  time  (Figure  3E).  The  goal 
of  this  deployment  was  to  capture  the  flow  velocity  while  the  USV  was  stationary  in 
order  to  compare  ADCP  measurements  in  high  discharge  environments  (ES)  and  low 
discharge  environments  (KR). 

C.  ELKHORN  SLOUGH,  CA,  AUGUST  2009 

Additional  velocity  profdes  obtained  from  bottom-mounted  upward-looking 
sensors  (Figure  5)  deployed  during  a  prior  experiment  that  lasted  9  days  in  August  2009, 
in  Elkhorn  Slough,  CA,  are  also  used  in  this  analysis  to  provide  true  stationary 
measurements  for  comparison  with  the  AV’s  ADCP  measurements.  Three  2  MHz  Nortek 
Aquadopp  ADCPs  and  one  1  MHz  Nortek  Acoustic  Wave  and  Current  (AWAC)  profiler 
were  deployed,  and  operated  continuously  with  a  sampling  rate  of  0.5  Hz  to  measure  the 
tidal  flow  in  the  water  column  with  a  bottom  blanking  distance  of  0.05  m,  0.35  m  bin 
size,  and  20  depth  bins.  These  stationary  ADCP  measurements  were  used  to  calculate 
averaging  windows  for  flows  ranging  from  0.40  to  1  m/s  for  comparison  with  the  AVs. 
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D.  BEAR  CUT  INLET,  FL,  JANUARY  2011 

The  UUV  was  utilized  in  January  2011,  in  Bear  Cut  inlet,  Miami,  FL,  referred  to 
as  BC.  BC  is  a  naturally  occurring  inlet  between  two  barrier  islands,  Virginia  Key  and 
Key  Biscayne  (Figures  3F  and  3G).  This  was  a  shallow  water  effort,  with  the  goal  of 
obtaining  velocity  profiles,  bathymetry,  and  water  quality  observations  in  an  environment 
experiencing  flows  greater  than  1  m/s.  Previously,  the  UUV  successfully  executed 
missions  in  flow  regimes  less  than  1  m/s.  The  UUV  was  used  in  “park”  mode  to  obtain 
stationary  velocity  measurements,  similar  to  the  station-keeping  USV  missions,  with  the 
UUV  facing  upstream  and  propelling  the  vehicle  as  necessary  to  stay  within  the  park 
radius  for  a  specified  amount  of  time.  At  BC,  the  UUV  parked  at  3  locations  for 
approximately  5  minutes  each,  collecting  ADCP  velocity  profiles.  The  ADCP  sampled  at 
1  Hz,  and  had  a  surface  blanking  distance  of  0.25  m,  bin  size  of  0.5  m,  with  30  depth 
bins. 


Location 

Flow 

Conditions 

(m/s) 

Depth 

(m) 

Platform 

Method 

Velocity 

Profiling 

Instrument 

Sampling 
Rate  (Hz) 

Number 

ofBins 

Bin 

Size 

(m) 

Blanking 

Distance 

(m) 

Conservative 

Mean  T* 

(min) 

Kootenai  River  - 
meandering  reach 
(August  2010) 

0.40 

8 

USV 

station- 

keeping 

EDI 

ADCP  (1.2 
MHz) 

0.5 

28 

0.25 

0.35 

2.9 

Kootenai  River  - 

braided  reach 

(August  2010) 

1.5 

3 

surface 

mim- 

catamaran 

fixed 

Nortek 
Aquadopp 
(2  MHz) 

1* 

35 

0.20 

0.05 

3.6 

Elkhom  Slough 

(Nov  ember  2010) 

0.40-1.00 

USV 

station- 

keeping 

EDI 

ADCP  (1.2 
MHz) 

0.5 

40 

0.25 

0.35 

2.8 

Elkhom  Slou^i 

(August  2009) 

0.40-1.00 

bottom  tri¬ 
pod 

fixed 

Nortek 
Aquadopp 
(2  MHz) 
and  Nortek 

AWAC 
(1  MHz) 

0.5 

20 

0.35 

0.05 

3. 2-3. 5 

Bear  Cut  inlet, 
Miami,  EL  (January 
2011) 

>1.00 

6 

UUV 

park 

mode 

Sontek 

ADCP 
(1.2  MHz) 

i* 

30 

0.5 

0.25 

■ 

Table  1.  Experiment  locations,  conditions,  instrument  platfonn  and  settings,  and 
conservative  mean  (mean  plus  one  standard  deviation)  averaging  window. 
Instruments  with  a  sample  rate  of  1  Hz,  which  were  averaged  to  a  new  0.5 
Hz  sample  rate  are  denoted  by  *. 
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Figure  3.  Experiment  locations:  (A)  Kootenai  River,  ID,  and  view  of  (B)  the 
meandering  reach  and  associated  USV  transects,  and  (C)  the  braided  reach  and  mini¬ 
catamaran  locations;  (D)  Monterey  Bay,  CA,  and  view  of  (E)  Elkhom  Slough  including 
the  USV  stationary  point  (green  dot)  and  the  stationary  ADCP  locations  (magenta  dots); 
(F)  Miami,  FL,  and  view  of  (G)  Bear  Cut  inlet  and  bathymetry  obtained  with  the  UUV 

and  the  UUV  park  locations 
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Figure  4.  Surface  mini-catamaran  with  downward-looking  ADCP. 


Figure  5.  Tri-pod  used  for  bottom-mounted  upward-looking  ADCPs 
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IV.  ANALYSIS  AND  RESULTS 


Examples  of  noisy  instantaneous  velocity  profiles  measured  by  the  station¬ 
keeping  USV  at  KR  and  by  the  UUV  in  “park”  mode  at  BC  are  shown  in  Figures  6 A  and 
6D.  In  order  to  reduce  the  statistical  noise,  time  averaging  is  required.  Profiles  of  velocity 
magnitude  averaged  over  the  duration  of  the  stationary  time  interval  (about  5  minutes) 
are  shown  in  Figures  6B  and  6E.  The  minimum  averaging  time  required  to  resolve  the 
instrument  and  environmental  noise  is  investigated  further.  By  averaging  measured 
velocity  profiles  over  increasing  sampling  times,  a  stable  estimate  of  the  mean  is 
eventually  reached  (Figures  6C  and  6F).  Previous  work  suggests  sampling  for  about  7  to 
1 1  minutes  at  a  fixed  location  based  on  visual  inspection  of  the  mean  velocity  as  a 
function  of  sampling  time  (Muste  et  al.  2004b).  However,  a  quantifiable  metric  to 
detennine  the  appropriate  averaging  time  is  desired. 

Two  statistical  methods  are  used  to  detennine  the  appropriate  averaging  window, 
T*,  such  that  a  stable  estimate  of  the  mean  flow  is  obtained:  1)  a  window-varying  time 
average,  and  2)  the  Kalman  Algorithm.  These  techniques  are  only  used  to  detennine  T* 
for  describing  the  appropriate  station-keeping  time,  and  are  not  meant  to  be  repeated  in 
the  field,  owing  to  the  need  of  prior  infonnation.  Both  methods  were  applied  to  a  10- 
minute  time-series  of  the  flow  magnitude  ( U  =  Vw2  +v2  )  and  direction  ( 0)  measured  by 
the  USV  with  a  dt= 2  s  in  the  KR  at  one  depth.  The  window-varying  time-averaging 
method  computes  the  mean  value  at  a  given  time  step  as  the  average  of  all  of  the  data  up 
to  that  time.  This  method  is  represented  by  the  following  equation: 

.«  population 

U(n)^=-  2>(l->")  ,  (2) 

n  , 

n=  1 

where  U  is  the  flow  magnitude,  n  is  the  number  of  observations  in  the  averaging  window, 
starting  at  n=  1  and  increasing  to  the  size  of  the  population.  This  equation  can  be  re¬ 
written  such  that  once  Uavg  is  calculated  for  a  given  averaging  window,  previous  data 
does  not  need  to  be  stored,  and  a  new  observation  simply  updates  Uavg  by 

U(n  + 1)  =—U(n  + 1),  +  —  U(n)  ,  (3) 

V  /  avg  fl  -\-\  ^  /observation  f%  ^  ' av^  v  7 
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where  U(n)avg  is  the  mean  from  the  previous  averaging  window  and  U(n+l)obsen,ation  is  the 
new  observation  at  n+ 1 .  Equation  3  shows  that  as  the  number  of  observations  increases, 
the  impact  of  a  new  observation  on  the  mean  decreases  by  1/(77  +1),  which  is  referred  to 
as  the  averaging  gain.  As  the  averaging  window  increases  the  velocities  asymptote  to  a 
constant  value,  which  is  the  mean  flow.  The  estimate  of  the  mean  U  and  0  becomes 
qualitatively  stable  for  averaging  windows  greater  than  60  s,  (Figures  7A  and  7C). 
However,  even  though  the  mean  appears  relatively  constant  after  60  s,  small  variations 
still  exist,  which  can  result  in  errors  when  describing  the  velocity  profile. 

A  statistical  metric  for  determining  when  the  time-averaged  signal  asymptotes  (an 
appropriate  T*)  is  required  such  that  the  observations  provide  a  stable  estimate  of  the 
mean  and  additional  observations  provide  minimal  new  information.  The  Kalman 
Algorithm  is  a  statistical  method  for  repeatedly  updating  the  estimate  of  the  mean  of  an 
evolving  system  from  a  sequence  of  “noisy”  measurements  by  processing  a  succession  of 
additional  measurements  (Kalman  1960).  The  Kalman  Algorithm  (KA)  is  defined  as: 

UN+l  =  K  x  (U observalion  —  UN^+  UN  ,  (4) 


where  Un  is  the  previous  estimate  of  U,  K  is  the  Kalman  gain,  Uobservation  is  the  new 
observed  value  at  ATI,  and  N  is  equal  to  the  number  of  independent  observations 
(degrees  of  freedom)  at  that  step.  K  is  defined  as: 


4/ 


K  = 


1 


V+fT2  l  +  N 

/  ^  w population 


as  N  — >  population 


(5) 


2  2 

where  a  n  is  the  variance  at  N,  and  a" population  is  the  system  variance,  which  is  assumed 
here  to  be  equal  to  the  variance  for  the  total  duration  of  10  minutes.  As  N  increases  to  the 
size  of  the  population,  cfN  =  crpopulation  and  K  =  l/(l  +  N).  When  K  is  small,  the  adjustment 


to  the  prediction  is  minimally  affected.  For  averaging  windows  less  than  60  s,  K  has  large 
initial  fluctuations,  but  then  exponentially  decays  for  a  larger  sample  size,  and  is  similar 
for  flow  magnitude  and  direction  (Figures  7B  and  7D). 

AT*  was  chosen  when  a  sufficient  number  of  independent  observations  were 
collected  to  obtain  a  K=0.01,  as  it  provides  a  conservative  threshold  for  describing  the 
mean  flow  magnitude  and  direction.  This  states  that  any  new  observation  only  improves 
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the  estimate  of  the  mean  by  1%  of  the  new  measurement  difference  (U0hservanon-UN).  As  K 
continues  to  decrease,  new  observations  provide  minimal  new  infonnation  to  the  system, 
even  if  the  new  observation  is  large.  A  larger  K  threshold  can  be  used,  resulting  in  a 
decrease  in  T*,  but  with  reduced  confidence  in  the  mean  velocity  estimate.  The  important 
aspect  of  this  approach  is  that  a  consistent  K  is  used  to  evaluate  ADCP  response.  Any 
reasonable  K  value  would  provide  satisfactory  results,  but  we  focus  on  T*  detennined 
when  A=0.01,  which  suggests  that  the  estimate  at  this  time  is  stable. 

Results  using  the  window-varying  time  average  method  (Equation  2)  and  the  KA 
method  (Equation  4)  are  comparable  (Figures  7A  and  7C).  The  largest  differences 
between  the  two  methods  occur  for  T*<60  s.  As  time  increases,  both  methods  are 
essentially  the  same  (Equation  3  and  Equation  4).  In  general,  KA  estimates  the 
asymptotic  mean  faster  than  just  averaging.  This  is  because  the  estimate  of  the  mean  by 
KA  (Equation  4)  uses  an  amount  proportional  to  the  difference  between  the  new 
observation  and  the  previous  estimate  of  the  mean  {U0bservati0n-UN),  whereas  the  averaging 
method  (Equation  3)  uses  the  entire  value  of  the  new  observation  (U observation),  which 
includes  more  noise  and  has  a  greater  impact  when  there  are  a  small  number  of 
observations.  Another  key  difference  between  the  two  methods  is  the  KA  retains  the 
system  variability  (Equation  5)  and  is  related  to  the  number  of  independent  observations, 
N,  while  the  window-varying  time  average  (Equation  3)  depends  solely  on  the  number  of 
observations,  n  (regardless  of  independence).  If  T*  was  estimated  when  the  averaging 
gain  equaled  0.01,  then  all  of  the  times  would  be  the  same,  as  it  is  strictly  dependent  on 
the  number  of  observations.  K  is  dependent  upon  the  population  variance  and  the 
instantaneous  variance,  which  will  vary  based  on  ADCP  manufacturer,  settings,  and 
environmental  conditions.  Differing  system  variance  for  varying  instrument  ADCP 
manufacturers  is  shown  in  Figure  8.  Regardless  of  manufacturer  the  variances  increase 
with  increasing  flow  magnitude.  The  KA  method  is  used  for  all  subsequent  analysis 
because  it  provides  a  quantifiable  metric  for  determining  the  appropriate  T*,  and  it 
accounts  for  the  variability  of  the  system  in  estimating  the  mean. 

The  criterion  for  the  decorrelation  time  is  when  the  autocorrelation,  R{f),  has 
decreased  from  the  initial  value,  R(z=  0  s),  by  a  factor  of  1/e.  The  decorrelation  time  of  U 
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measured  by  both  the  USV  in  KR,  with  dt= 2  s,  and  an  Aquadopp  in  KR,  with  dt=  1  s, 
occurs  in  one  dt  for  each  instrument  (Figure  9)  and  is  consistent  with  Muste  et  al. 
(2004b).  For  comparison,  measurements  taken  by  an  Acoustic  Doppler  Velocimeter 
(ADV),  deployed  simultaneously  in  the  braided  reach  of  KR,  were  examined  because 
ADVs  can  sample  faster  (dt=\l2>2  s).  The  UAdv  signal  was  low-pass  (LP)  fdtered  with  a 
frequency  cut-off  of  1  Hz,  which  was  the  onset  of  the  noise  floor.  Uadv,lp  decorrelates  in 
less  than  1  s,  suggesting  that  the  environmental  fluctuations  are  temporally  short  (Figure 
9).  ADCPs  intrinsically  sample  at  a  faster  rate  (<1  s),  but  record  an  average  of  these  faster 
observations  at  a  slower  rate  (>1  s),  which  is  longer  than  the  environmental  decorrelation 
time  as  found  by  the  ADV.  Therefore,  a  decorrelation  time  is  limited  by  the  ADCP’s 
slower  sample  “recorded”  rate  (>1  s).  ADCPs  with  bottom-tracking  have  an  even  slower 
sample  rate  (>  2  s)  because  the  ADCP  is  sampling  both  water  profile  and  bottom  track 
estimates  independently.  The  Aquadopps  and  AWAC  do  not  have  bottom-track 
capabilities,  allowing  for  a  faster  “recorded”  sample  rate.  The  focus  of  the  manuscript  is 
for  ADCPs  with  bottom-tracking,  as  these  are  the  systems  commonly  found  on  AVs. 
Therefore,  the  re-sampling  of  all  instruments  to  a  dt= 2  s  eliminates  the  influence  of  dt 
when  evaluating  the  effects  of  depth,  environmental  conditions,  platform,  and  instrument 
manufacturer  on  fi  and  provides  a  conservative  estimate  for  the  number  of  independent 
observations.  Note  that  the  relative  flow  magnitude  and  direction  of  the  environment  was 
obtained  in  approximately  60  s,  when  dt=  2  s  (Figures  7A  and  1C).  This  corresponds  to  30 
observations,  which  is  considered  a  large  sample  size  and  is  statistically  a  sufficient 
number  to  accurately  describe  the  mean  (Davis  2002).  However,  Muste  et  al.  (2004b) 
showed  that  the  time  requirement  to  resolve  the  velocity  mean  exceeds  this  statistical 
parameterization.  Measurements  acquired  in  differing  environments  with  several 
instruments  and  platforms  are  evaluated  to  determine  a  statistically  accurate  averaging 
time. 

The  relationship  between  T*  computed  for  the  flow  magnitude  (U)  and  for  the 
flow  direction  (9)  using  depth- averaged  measurements  from  the  station-keeping  USV  in 
KR  is  poor  (R=0.18)  suggesting  there  is  no  dependent  relationship  (Figure  10). 
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However,  both  measures  indicate  T*  of  approximately  2.5  minutes.  In  other  words,  the  T* 
for  the  cross-  and  along-channel  velocities  are  independent,  but  similar. 

Velocity  measurements  from  all  of  the  experiments  were  analyzed  to  detennine 
how  T*  varies  by  depth,  environment,  platfonn,  and  instrument  manufacturer.  The  mean 
T*  values,  as  a  function  of  depth,  calculated  from  stationary  measurements  by  the 
differing  platfonns  and  ADCPs  during  each  experiment  are  shown  in  Figure  1 1 .  Based  on 
a  95%  confidence  level,  T*  is  only  statistically  uniform  over  the  vertical  for  the  USV 
measurements  in  ES  and  for  the  mid-water  column  measurements  by  the  Aquadopps  and 
AWAC  in  ES.  However,  the  T*  values  are  relatively  the  same  throughout  the  water 
column  for  all  of  the  measurements,  with  the  exception  of  the  UUV,  regardless  of  the 
platform,  ADCP  manufacturer,  and  environmental  conditions.  Only  three  vertical  profiles 
were  measured  with  the  UUV,  therefore  there  were  not  enough  data  to  conclude  that  the 
mean  vertical  profile  was  statistically  uniform. 

Examining  the  results  from  the  A  Vs  in  differing  environments  suggests  T*  is 
similar  regardless  of  flow  conditions  for  each  platfonn,  with  relatively  small  differences 
in  T*  between  AVs.  In  a  low  flow  environment  (0.4  m/s)  in  the  meandering  reach  of  the 
KR,  T*  for  the  station-keeping  USV  varied  throughout  the  river  resulting  in  a  mean  T*  of 
approximately  2.6  minutes  (Figure  11  A).  The  standard  deviation  is  0.3  minutes  resulting 
in  a  conservative  T*  limit  of  2.9  minutes.  Comparatively,  in  a  higher  flow  environment 
(0.4  -  1  m/s)  in  ES,  the  mean  and  standard  deviation  for  T*  for  the  station-keeping  USV 
was  similar  at  2.5  and  0.3  minutes,  respectively  (Figure  11B).  The  results  are  the  same 
for  the  USV  in  the  varying  flow  conditions.  Owing  to  these  values,  it  is  recommended 
that  the  USV  station-keep  at  the  same  location  for  a  minimum  of  3-4  minutes  to  ensure 
that  a  stable  estimate  of  the  mean  is  obtained.  In  the  fastest  flows  (>1  m/s)  in  BC,  the 
mean  and  standard  deviation  of  T*  measured  by  the  stationary  UUV  (Fig.  1 1 C)  were 
similar  to  the  station-keeping  USV  results,  at  3.2  and  0.1  minutes,  respectively. 

Although  the  AVs  were  considered  stationary  when  performing  station-keeping 
and  “park”  missions,  there  was  still  some  vehicle  movement.  The  speed  of  the  USV  while 
station-keeping  was  negligible  (<  0.01  m/s)  and  the  UUV  moved  an  average  of  0.17  m/s 
while  in  park  mode.  The  speed  of  the  AVs  while  maintaining  a  stationary  position  was 
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determined  based  on  bottom-track  velocities,  which  have  less  than  a  6%  error  compared 
to  kinematic  GPS  measurements  (Fong  and  Monismith  2004).  To  further  evaluate 
potential  errors  that  developed  because  of  the  moving  platform,  T*  for  the  bottom- 
mounted  stationary  ADCPs  deployed  in  ES  were  estimated  (Figures  1  IE  and  ,1  IF).  Note 
that  these  stationary  ADCP  measurements  were  not  collected  at  the  same  time  as  the 
USV  mission  in  ES.  The  conservative  T*  (mean  plus  one  standard  deviation)  for  the 
bottom-mounted  Aquadopps  and  AWAC  were  3.2  and  3.5  minutes,  respectively.  These 
T*  values  are  comparable  to  the  station-keeping  USV,  even  though  the  ADCP 
manufacturers  are  different.  These  results  indicate  that  the  station-keeping  A  Vs  are 
stationary  enough  to  obtain  stable  estimates  of  the  mean  flow  in  similar  time  frames  as 
true  stationary  measurements.  Additionally,  the  downward-looking  Aquadopps  mounted 
on  stationary  surface  mini-catamarans  in  KR  (Figure  11D)  resulted  in  a  conservative 
mean  T*  of  3.6  minutes,  which  is  in  the  range  of  the  other  measurements. 

The  depth-averaged  T*  for  the  various  platforms  in  different  flow  environments 
are  shown  in  Figure  12.  The  mean  T*  was  found  to  be  about  3  minutes  with  a  maximum 
T*  of  less  than  4  minutes,  with  a  weak  relationship  with  the  depth-averaged  velocity.  T* 
for  the  USV  measurements  varied  by  about  1  minute  over  a  relatively  small  velocity 
range,  whereas  T*  for  the  Aquadopp  in  KR  was  constantly  about  3.2  minutes  for  a  larger 
velocity  range  (1.5  m/s).  However,  the  Aquadopps  deployed  in  ES  had  T*  values  ranging 
from  0.5  to  3.5  minutes  for  a  similar  large  velocity  range.  Therefore,  there  is  no 
significant  dependence  of  T*  on  flow  condition,  platform,  or  ADCP  manufacturer,  and  a 
conservative  value  of  4  minutes  will  provide  a  stable  estimate  in  any  scenario. 
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Figure  6.  ADCP  data  collected  by  the  station-keeping  USV  in  KR  (top  row)  and  the 
UUV  in  “park”  mode  at  BC  (bottom  row):  (A,D)  raw  flow  magnitude  as  a  function  of 
time  and  depth,  (B,E)  time-averaged  vertical  flow  magnitude  profdes,  and  (C,F) 
window-varying  time  averaged  flow  magnitude  as  a  function  of  averaging  window  and 
depth.  Colorbars  represent  flow  velocity  magnitude. 
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Figure  7.  Values  of  (A)  flow  velocity  magnitude  estimated  by  the  window-varying 
time  average  method  (black)  and  by  KA  (blue),  (B)  velocity  magnitude  Kalman  gain,  K, 
(C)  flow  direction  estimated  by  the  window-varying  time  average  method  (black)  and  by 
KA  (blue),  and  (D)  direction  Kalman  gain,  K,  as  a  function  of  averaging  window,  T*,  for 
one  ADCP  bin  from  the  station-keeping  USV  in  KR.  The  red  circles  indicate  when 
K<0.01  and  the  estimate  of  the  mean  is  considered  stable. 
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Figure  8.  Depth- averaged  flow  velocity  magnitude  variance  versus  depth- averaged 

mean  flow  velocity  magnitude  measured  by  ADCPs  mounted  on  various  stationary 

platforms  in  differing  environments. 


Figure  9.  Autocorrelation  functions  as  a  function  of  time  lag,  t,  for  U  measured  by: 

the  station-keeping  USV  in  KR  with  dt=  2  s  (black),  an  Aquadopp  in  KR  with  dt=  1  s 
(blue),  and  the  ADV  in  KR  with  dt=  1/32  s  (green),  where  the  red  dashed  line  represents 

the  decorrelation  criteria  (1/e). 
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Figure  10.  Correlation  between  T*  calculated  using  flow  velocity  magnitude  ( U)  and 
T*  calculated  using  flow  direction  (0).  The  black  solid  line  represents  an  idealized  linear 
regression  line  and  the  red  solid  line  represents  the  actual  liner  regression  line  for  the 

comparison  (R2=0.18,  m=0.43). 
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Figure  11.  Averaging  window,  T*,  as  a  function  of  depth,  detennined  for  (A)  the 
station-keeping  USV  at  KR,  (B)  the  station-keeping  USV  at  ES,  (C)  the  UUV  in  “park” 
mode  at  BC,  (D)  stationary  downward-looking  Aquadopps  at  KR,  and  for  the  fixed 
bottom-mounted  (E)  Aquadopps  and  (F)  AWAC  in  ES.  Note  that  all  ADCPs  were 
stationary.  Solid  lines  depict  the  mean  averaging  window,  dashed  lines  depict  one 
standard  deviation,  and  dotted  lines  depict  the  95%  confidence  level. 
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Figure  12.  Depth- averaged  flow  velocity  magnitude  versus  depth-averaged  averaging 
window,  T*,  measured  by  ADCPs  mounted  on  various  stationary  platforms  in  differing 

environments. 
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V.  SUMMARY 


As  with  any  observational  endeavor,  an  acceptable  prior  level  of  statistical 
confidence  is  required.  The  effectiveness  of  using  ADCPs  on  AVs  at  fixed  locations  to 
obtain  accurate  measurements  of  the  mean  flow  conditions  in  various  rivers  and  estuaries 
has  been  investigated.  Using  the  Kalman  Algorithm  with  a  Kalman  gain  of  K=0.01 
provided  a  useful  method  for  estimating  a  statistically  relevant  T*  for  obtaining  high- 
quality  ADCP  horizontal  velocity  profiles.  The  mean  T*  for  various  ADCPs  and 
platfonns  was  found  to  be  3  minutes  with  a  maximum  T*  of  less  than  4  minutes,  with 
minimal  dependence  upon  instrument  type  and  the  environmental  conditions.  Surprising 
to  the  authors,  all  ADCPs  and  their  platforms  responded  similarly.  The  conservative  T*  is 
two  times  smaller  than  those  found  in  previous  studies  examining  moving  and  fixed 
vessel  ADCP  measurements  (Muste  et  al.  2004b,  Szupiany  et  al.  2007).  The  instrument 
sampling  rate,  dt,  can  influence  T*,  with  T*  decreasing  with  decreasing  dt.  T*  is  not 
dependent  on  depth  or  flow  velocity,  and  there  is  no  correlation  (R'=0.18)  between  using 
flow  velocity  magnitude  and  direction.  A  conservative  time  of  4  minutes  of  observations 
at  a  specific  point  should  be  acquired  to  obtain  an  accurate  estimate  of  the  mean  velocity 
profile.  If  ADCP  observations  are  collected  with  a  moving  vessel,  the  velocities  within  the 
internal  travelled  in  4  minutes  should  be  spatially  stationary.  The  statistical  techniques 
implemented  in  this  work  provided  a  robust  means  of  detennining  the  4-minute  averaging 
time  window  for  these  platforms,  ADCPs,  and  riverine/estuarine  environments.  The 
techniques  do  not  need  to  be  repeated,  unless  these  conditions  significantly  change. 

The  UUV  was  able  to  perform  surface  station-keeping  (“parks”),  however,  the 
UUV  works  best  when  operating  below  the  surface  and  at  a  speed  greater  than  1  m/s.  At 
these  speeds,  the  UUV  will  have  traveled  approximately  240  m  in  the  necessary 
averaging  window.  If  flow  along  this  transect  can  be  considered  temporally  and  spatially 
stationary  (homogeneous),  an  accurate  description  of  the  profile  can  still  be  resolved. 

The  USV  can  precisely  execute  repeated  station-keeping  positions  in  varying 
environments.  An  example  of  the  USV’s  spatial  velocity  measurement  capability  is 
shown  in  Figure  13.  In  planning  the  mission,  a  station-keeping  time  of  10  minutes  was 
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used  to  resolve  the  mean  velocity  profiles  (Muste  et  al.  2004b).  The  USV  successfully 
and  efficiently  mapped  the  velocity  structure  (Figure  13)  in  the  meandering  reach  of  the 
KR  in  2  days.  With  the  new  knowledge  that  a  conservative  T*  of  4  minutes  is  necessary, 
the  number  of  observational  locations  could  have  doubled,  resulting  in  an  even  better 
description  of  the  river  flow  field.  Consequently,  a  better  understanding  of  a  river’s  flow 
field  will  provide  the  U.S.  Navy  with  the  ability  for  better  mission  planning.  The 
capabilities  of  the  AVs  to  map  the  3D  flow  conditions  of  an  environment  accurately  and 
efficiently,  in  time  and  space,  far  exceed  the  capabilities  of  human-controlled  vessels. 
Because  of  this  fact,  AVs  are  of  particular  interest  to  the  U.S.  Navy  because  they  increase 
the  safety  of  its  personnel  by  reducing  the  human  footprint  that  was  necessary  to  collect 
this  type  of  environmental  data  in  the  past.  AVs  provide  scientists  an  indispensable  tool 
to  effectively  study  various  flow  environments.  Additionally,  AVs  provide  the  U.S.  Navy 
with  environmental  intelligence  that  will  enhance  overall  mission  success  and  a  means  to 
protect  its  personnel  from  prolonged  exposure  to  any  threats,  human  or  environmental, 
that  would  cause  them  harm  . 


28 


latitude 


longitude 

Figure  13.  Spatial  map  of  the  velocity  field  in  the  meandering  reach  of  the  KR 
measured  by  the  station-keeping  USV. 
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